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The intergalactic magnetic field (IGMF) in cosmic voids can be indirectly probed through its effect on electro- 
magnetic cascades initiated by a source of TeV gamma rays, such as blazars, a subclass of active galactic nuclei. 
Blazars that are sufficiently luminous at TeV energies, "extreme TeV blazars" , can produce detectable levels of 
secondary radiation from inverse Compton scattering of the electrons in the cascade, provided that the IGMF 
is not too large. We reveiw recent work in the literature which utilizes this idea to derive constraints on the 
IGMF for three TeV-detected blazars-lES 0229+200, 1ES 1218+304, and RGB J0710+591, and we also inves- 
tigate four other hard-spectrum TeV blazars in the same framework. Through a recently developed detailed 
3D particle tracking Monte Carlo simulation code, incorporating all major effects of QED and cosmological 
expansion, we research effects of major uncertainties such as the spectral properties of the source, uncertainty 
in the intensity of the UV - far IR extragalactic background light (EBL), under-sampled Very High Energy 
(VHE; energy > 100 GeV) coverage, past history of gamma-ray emission, source vs. observer geometry, and jet 
AGN Doppler factor. The implications of these effects on the recently reported lower limits of the IGMF are 
thoroughly examined to conclude that presently available data are compatible with a zero IGMF hypothesis. 



1. Introduction 

It is well known that astrophysical magnetic fields 
are ubiquitous in galaxies and galaxy clusters on the 
order of a micro Gauss 
[2001 



(see e.g. 



Grasso and Rubinstein 



Widrow 2002 ). Furthermore, there is increas- 
ing theoretical evidence based on numerical simula- 
tions of structure formation that nano Gauss order 
fields permeate filaments of the large scale structure 
(see e.g. |Ryu et al. 2008| ). However, at present there 
has been no detection of the intergalactic magnetic 
field (IGMF) presumed to exist in the Cosmic Void 
regions of the large scale structure. The detection of 
the IGMF could provide important insights for solv- 
ing outstanding problems of its origin and role in both 
the cosmology and astrophysics of structure formation 
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Until recently, only upper limits on Bjqmf have 
been established, on the order of Bjqmf ~ 10 -9 
Gauss, at a correletion length of 1 Mpc, and this limit 
weakens as the correlation length (Ac ) decreases as 



A, 



-1/2 



Neronov and Semikoz 



2009 



However, a 



more sensitive technique has emerged during the past 
few years which may become a tool for the measure- 
ment of IGMF characteristics. This technique relies 
on observations of blazars, in the energy range from 
100 McV to greater than 10 TeV and is described 



by several authors Neronov and Semikoz 2006 Eu 



ngwanichayapant and Aharonian 2009, Dolag et al 



2009[|Elyiv et al.|2009| . Briefly, TeV-scale gamma rays 
from the blazar interact with the UV to far-IR extra- 
galactic background light (EBL), producing electron- 
positron pairs, which then undergo inverse Compton 
(IC) scattering on the CMB photons, producing sec- 
ondary gamma rays of a lower energy than the pri- 
mary. Because the pairs' trajectories depend on the 
interactions with the magnetic field, comparing the 
GeV-scale secondary gamma ray radiation with exist- 



ing data can be used to characterize properties of the 
IGMF. 



In the past few years, a number of studies Neronov 


and Vovk |2010| Tavecchio et al.| 20101 


Dermer et al. 


20111 |Huan et al. 2011 Dolag et al. 


2011 Taylor 


et al. 2011 have demonstrated that a 


ower limit on 



the IGMF strength can in principle be derived by 
requiring that the secondary gamma ray GeV radi- 
ation which would be produced at a given Bjgmf 
not exceed the measured value. A consensus has be- 
gun to emerge that a minimum value of Bjgmf ~ 
10 -17 — 10~ 18 Gauss is required to explain the obser- 
vations. With a newly developed 3-dimensional Monte 
Carlo code combined with a detailed and general sta- 
tistical analysis method, the purpose of these proceed- 
ings is to show that the current observational data is 
consistent with the Bjqmf — hypothesis, when the 
systematic and statistical uncertainties are accounted 
for. 



2. Monte Carlo Simulations 

In order to explore in detail the potentially observ- 
able effects of cascading in the voids of the large scale 
structure, we have developed a fully 3-dimensional 
Monte Carlo code, which propagates individual parti- 
cles of the cascade in a cosmologically expanding Uni- 
verse and accounts for the QED interactions with the 
EBL and CMB without simplifications. 

The spectral energy density of the EBL from the 
UV to the far-IR is not precisely known, and therefore 
we choose to model it with specific SEDs close to the 
energy density of the recent EBL model of |Dommguez| 
The magnetic field is described as a 



et al. 2011 



system of cubic cells with an edge length set to 1 Mpc, 
with equal field amplitudes in each cell, but randomly 
oriented in direction. 

The gamma ray source model employed is based on 
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leading theoretical speculations about the nature of 
the TeV blazar source (see e.g. Urry and Padovani 



1995 ). In the reference frame of the blazar jet, GeV 



TeV photons are distributed isotropically with a bro- 
ken power law spectrum of index a for the higher en- 
ergy component and 7 for the lower energy compo- 
nent. Once this distribution is boosted into the ref- 
erence frame of the host galaxy with a doppler boost 

— 1/2 

factor of T = (l — /3 2 ) , the observed differential 
flux energy density (dFED) can be parameterized as 



dF 
de 



F n 6 2 




-7+1 



exp 



Q + l 



exp 



(?) 
(?) 



< 1 



> 1 



where e is the photon energy in the reference frame 
of the host galaxy, S = [r(l — /3cos0 v )]~ > $v is the 
viewing angle from the blazar jet axis to the line of 
sight of the observer, F is a flux normalization factor, 
e c is the exponential cutoff energy, and cb is the spec- 
tral break energy. This six parameter 7-ray source 
spectrum is given at the redshift of the host galaxy 
and is necessary and sufficient to satisfy observational 
data of TeV blazars in both the HE (Fermi-LAT) and 
VHE (IACT) regimes, and it is more general than 
previous studies which only considered single intrinsic 
power law models for the primary gamma ray emis- 
sion, over the energy range spanning more than five 
orders of magnitude from ss 100 MeV to > 10 TeV. 



3. Testing the B 



IGMF 



hypothesis 



The most comprehensive work yet on constraining 
the IGMF was done by |Taylor et aT] [2011] (here- 



after referred to as TVN11), in which three sources 
were found (RGB J0710+591, 1ES 1218+304, and 
1ES 0229+200) whose combined quasi-simultaneous 
IACT and Fermi-LAT dFED are not compatible with 
the Bjgmf — hypothesis at greater than 95% con- 
fidence level. It was suggested that Bjgmf > 10~ 17 
G is required to explain observations. 

With the newly developed simulation and analy- 
sis code, we re-examined the conclusion of a lower 
limit on Bjgmf- The spectral model of Eq [T] was 
fit to the same IACT data as used in TVN11, but 
with about one more year of Fermi data. A large 
range of parameter space of the model was scanned 
over (a, e c , F Q , 7, £b), while several parameters were 
left fixed, including 1) a fixed source-observer geom- 
etry of T = 10 and V = 0°, 2) the EBL mod el with 
an SED very close to Dommguez et al. [201 1] , 3) the 
duty cycle of the TeV emission, and 4) assuming no 
structures with an enhanced magnitude of magnetic 
field (for ex: galaxy clusters, filaments, etc.) between 
the source and observer. For each model, a x 2 fit to 
the combined IACT /Fermi data was performed, and 



Arlen et al 



source are s 



also the corresponding confidence level for rejecting 
the Bjgmf = hypothesis. (For more details on this 
procedure, as well as the chosen data sets, please see 
2012| ). The best fitting models for each 
lown in Fig[l] The first two sources, RGB 
J0710+591 and 1ES 1218+304 (Fig[l]a and b) were 
found to be incompatible with the Bjgmf — hy- 
pothesis at < 90 %, which does not corroborate the 
conclusion of TVN11 that the Bjgmf = hypothe- 
sis is rejected for these sources. It was found that a 
combination three factors were the primary reasons 
for this discrepancy, which are the following: 1) the 
updated Pass 7 Fermi-LAT data used in the present 
work, 2) the more general description of the broken- 
power law source model in the analysis of the data, 
and 3) the more robust statistical analysis which used 
information about the cascade flux to define the spec- 
tral index in each bin as an input to the Fermi likeli- 
hood spectral analysis. 

It is confirmed that the data for 1ES 0229+200 
are incompatible with the Bjgmf = hypothesis at 
> 99% confidence level, with the standard assump- 
tions of 1) - 4) (listed in the previous paragraph). 
Since only this single source was found to be in con- 
flict with a non-zero IGMF, each of these standard 
assumptions were varied to within their uncertainty 
limits to research all the possible ways in which the 
non-zero IGMF requirement can fail. For example, the 
duty cycle assumption was varied, because the 3 year 
time-averaged Fermi-LAT spectrum is combined in 
the spectral model fits with the VHE spectrum which 
is averaged over a much shorter time period. Thus, 
it is necessary to make an assumption regarding the 
stability of the pointed VHE observations (duty cy- 
cle). To investigate the effect of a reduced duty cy- 
cle, the spectrum of this source was modified at the 
highest 5 energy points to half an order of magni- 
tude of their reported values. It was found that the 
VHE-HE data set combined in this way does not rule 
out the Bjgmf = hypothesis to more than 95 % 
confidence level. Therefore, if the true 3-year time- 
averaged spectrum is over-estimated by about 50 %, 
then compatibility with the Bjgmf = hypothesis 
would be achieved. A similar compatibility could be 
achieved if the location of the source within a filament 
is separated by at least a few hundred Mpc from the 
nearest Cosmic Void. The enhanced magnetic field 
within a filament could potentially isotropize the elec- 
trons before IC scattering occurs along the line of sight 
to the obsvcrvcr and could also explain the observa- 
tions without invoking a non-zero IGMF. 

Alternatively, the energy density of the EBL in the 
near IR band can be reduced to the lower limits de- 
rived from galaxy counts data to avoid incompatibility 
with the zero IGMF assumption. This latter case of 
EBL uncertainty has also been investigated recently 



Vovk et al. 112012 to conclude that this data provides 



joint constraints on the IGMF and EBL in which the 
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Figure 1: Best fit models for (a) (left) source RGB J0710+591 where a = 1.8, tc = 3.16 TeV, which is inconsistent 
with the Big m f = hypothesis at < 90% (b) (middle) source 1ES 1218+304 where a — 1.8, and ec = 3.16 which is 
inconsistent with the Bigmf = hypothesis at < 80% (c) (right) source 1ES 0229+200 where a = 1.3, ec = 1 TeV, 
which is inconsistent with the Bigmf = hypothesis at > 99%. 



2.5 



1.5 



a 1 



0.5 




- prompt+cascade 
■ prompt 

cascade 

- IACT data 

- Fermi LAT data 



Energy Cut [GeV] 




10 10' 
Energy [GeV] 



Figure 2: (a) (left) Confidence level for the simulated source models for 1ES 0229+200 obtained under the assumption 
of a low EBL Model close to the lower limits due to galaxy counts. The black line terminated with crosses represents 
the range of models analyzed in |Vovk et af][2012| . (b) (right) dFED of 1ES 0229+200 for the best fitting model at 
a = 1.6, e c = 3.16 TeV. 



zero field is in the exclusion region. Figure[2]illustrates 
our finding that this conclusion is the result of a lim- 
ited exploration of the parameter space of uncertain- 
ties in the intrinsic properties of the source. There- 
fore, for EBL models which are very low in the ncar-IR 
energy density but still allowable by the currently es- 
tablished lower limits, the data of 1ES 0229+200 can 
be reconciled with a zero magnetic field (Bigmf = 0). 



4. Discussion 

We have investigated the HE-VHE energy spec- 
trum of several extreme TeV blazars for which radi- 
ation in the HE band may be dominated by the sec- 
ondary photons produced through cascading in the 
voids along the line of sight. These sources are char- 
acterized by their observed hard spectra in the VHE 
band accompanied by redshifts of order ss 0.1 suggest- 
ing a very large energy output into pair production 
and subsequent cascading. For these sources, obser- 



vations in the HE band can therefore limit the flux 
of secondary photons and establish a lower bound on 
the IGMF. This strategy has been utilized in several 



publications Neronov and Vovk|2010 


Tavecchio et al. 


2010[ |Dermer et al. 2011 Huan et al. 2011 Dolag 


et al. 2011 


Taylor et al. 


2011 to sug 


;gest Bigmf < 



In contrast to these studies, we systematically inves- 
tigated effects of a wide range of uncertainties using 
detailed 3D Monte Carlo simulations to conclude that 
the Bigmf — hypothesis remains compatible with 
current observations. 
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